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We investigate the decoherence process of three qubit system obtained by manipulating the state 
of a trapped two-level ion coupled to an optical cavity. Interaction of the ion with a resonant 
laser and the cavity field tuned to red sideband of ionic vibrational motion, generates tripartite 
entanglement of internal state of the ion, vibrational state of ionic center of mass motion and the 
cavity field state. Non-dissipative decoherence occurs due to entanglement of the system with the 
environment, modeled as a set of noninteracting harmonic oscillators. Analytic expressions for the 
state operator of tripartite composite system, the probability of generating maximally entangled 
GHZ state, and the population inversion have been obtained. Coupling to environment results in 
exponential decay of off diagonal matrix elements of the state operator with time as well as a phase 
decoherence of the component states. 

Numerical calculations to examine the time evolution of GHZ state generation probability and 
population inversion for different values of system environment coupling strengths are performed. 
Using negativity as an entanglement measure and linear entropy as a measure of mixedness, the 
entanglement dynamics of the tripartite system in the presence of decoherence sources has been 
analyzed. The maximum tripartite entanglement is found to decrease with increase in the strength 
of system-environment coupling. The negativity as well as the linear entropy as entanglement 
measures give qualitatively similar results, uniquely identifying maximally entangled and separable 
states of the system. For large values of system-environment coupling strength, the mixed states of 
the composite system lying at the boundary of entangled-separable region are reached. For these 
states the negativity and linear entropy show distinctly different behaviour. This can be understood 
by noting that whereas the negativity measures entanglement generating quantum correlations, the 
linear entropy measures all correlations that reduce the purity of the state. 

Keywords: Trapped ions, Cavity QED, Tripartite entanglement, GHZ state generation, Non- 
Dissipative Decoherence, Negativity, Linear entropy, Mixed states 



INTRODUCTION 

Controlled manipulation of quantum states and implementation of quantum logic gates are essential elements 
of a quantum computer. Quantum computation relies on nonclassical properties of qubits such as entanglement. 
Entanglement dynamics is of utmost importance to quantum communication dense coding 0] and quantum 
teleportation Q, as well. Cold ions in a linear trap offer a promising physical system to implement quantum 
computation, as each ion allows two qubit state manipulation. One of the experimental efforts to realize multipartite 
entanglement involves trapping two level ion in an optical cavity. Coupling of trapped ion to quantized field inside 
an optical cavity has been successfully achieved More than one trapped ions can be used to construct a more 
complex multiple function system in which quantum states of several trapped ions entangled with cavity photons are 
manipulated in a controlled manner. A string of trapped ions in a cavity may be used for implementing quantum 
gates . In quantum networks involving cavity QED setups Q , quantum information is stored and processed by ions 
trapped in a cavity. De-coherence of quantum systems is one of the major hurdles to implement quantum computers. 
Suppression of quantum coherence occurs due to interactions of the quantum system with environment resulting 
in random and unknown perturbations of system Hamiltonian. In a previous article ||, we proposed a scheme to 
generate three qubit maximally entangled GHZ state, using trapped ion interacting with a resonant external laser 
and sideband tuned single mode of a cavity field. In this article, we investigate the decoherence process of three qubit 
system obtained by manipulating the state of a two level trapped ion coupled to an optical cavity. In the context of 
a qubit based on a single Ca + ion, Schmidt-Kaler et al |9( have found magnetic field fluctuations and laser frequency 
fluctuations to be the major decoherence sources. Decoherence of motional quantum states of a trapped atom coupled 
to engineered phase and amplitude reservoirs has been measured [ToL . The trap frequency was changed by applying 
potentials to trap electrodes to simulate a phase reservoir. In general the environment or heat bath is modeled as a 
system of noninteraction boson modes 01 • It is known that for harmonic systems a heat bath is effectively equivalent 
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to an external uncolrrelated random force acting on the quantum system. As such the decoherence effects due to 
magnetic field fluctuations, laser frequency fluctuations, and potentials applied to elctrodes can be modeled as a 
system of noninteracting boson modes. Other important decoherence sources are cavity losses and spontaneous decay. 
The case where spontaneous decay effects and cavity losses become important has been examined by Fidio et al. [13| . 
In the present work, these two effects are assumed to be negligible. 

As the laser-ion interaction times involved in experiments with ions trapped in a cavity are in the micro seconds 
range, Markovian approximation is likely to yield undesirable results. We consider here only adiabatic decoherence of 
the system with no energy exchange between the quantum system and environment. A detailed discussion on adiabatic 
decoherence has been given by Mozyrsky and Privman in ref . 0] . Starting from an initial state in which the system 
and environment are in a separable state, the state operator of the ion-cavity system is obtained by tracing over the 
environment degrees of freedom. Coherent states are used to evaluate the trace [L^ . The decoherence effects result in a 
state operator with off diagonal matrix elements decaying with time, in the eigen basis of the interaction Hamiltonian. 
The state operator is used to examine the decoherence effects on the probability of generating maximally entangled 
tripartite GHZ state and population inversion of two level ion. As the coherences decrease the system tends to a 
purely mixed state of energy eigenvectors. Decoherence of tripartite entanglement is analyzed by applying the Peres- 
Horodecki condition for separability 0] to bipartite decompositions of the composite system. A state is separable if 
the partial transpose of it's density matrix with respect to a subsystem is positive semidefinite. An entangled state 
violates the positive partial transpose (PPT) separability criterion. Negativity is an entanglement measure based 
on (PPT) separability criterion. We use Negativity to characterize the entanglement and linear entropy to measure 
the mixedness of three subsystems. 



UNITARY EVOLUTION OF ISOLATED SYSTEM 



Consider a single two-level ion trapped in a high finesse optical cavity. The frequency of radio frequency (Paul) 
trap along x axis is v and the excitation energy of the ion is Two®. The ion is radiated by an external resonant laser of 
frequency ljl an d interacts with the single mode cavity field of frequency to c tuned to red sideband of ionic vibrational 
motion. Interaction with external laser field as well as the cavity field generates entanglement of internal states of the 
ion, vibrational states of ionic center of mass, and the cavity field number state. The free Hamiltonian of the system 
is 

H os — \$a + -J + huj^b + —^-^z , (1) 

cr (a) and b'(b) being the creation(destruction) operators for vibrational phonon and cavity field photon respectively. 
With the ion placed close to the node of cavity field standing wave, the interaction Hamiltonian is given by 

Hi = M7[<7+ exp \irn{a + a') — iu>i,t\ + <r_ exp [ir]i(a + a') + iwi,t\ 

+hg (a+b + a-tf) sin [r? c (a + a f )] , (2) 

where f2 , g, til, f] c are the Rabi frequency, the ion-cavity coupling strength, ion-laser Lamb Dicke parameter, and 
ion-cavity field Lamb Dicke parameter, respectively. We work in the Lamb Dicke regime that is t]l -C 1 , r\ c <C 1 . In 
the interaction picture determined by the transformation U — exp(iHost/h) and rotating wave approximation, for 
u>l = wo and tu c = luq — v, the interaction Hamiltonian reduces to 

Ha = Hfl[a + + <r_] + figrj c a+ba + aJb* a) . (3) 

The unitary evolution of the isolated system in interaction picture is given by 

ih d -Ml=Hii*i(t), 
at 



where ^>i(t) = exp(iH st/h)^s(t)- 
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Basis truncation 



We expand Hji in the space of basis vectors \i,k,l) as 

II n = > (i',k',l'\H H \i,k,l)\i',k',l') (i,k,l\ , 



E 

i,k,l,i' ,fc' 



(4) 



where i = (gore) and fc, 1(0, 1, .., oo), denote the state of ionic vibrational motion and the cavity field num- 
ber state, respectively. We isolate a four dimensional vector space containing the computational basis vectors 
\g, m — 1, n — 1) , |e, to — 1, n — 1) , |<?, to, n) and |e, to, n) for a given choice of to, n values. The matrix represent- 
ing the interaction Hamiltonian acting in four dimensional space is 



His 






Ml 








nn 





hgr/c^/mn 













nn 








Tift 






(5) 



The unitary transformation that diagonalizes His is easily obtained and yields the eigenvectors satisfying His$ p 
EpQp, (p = 1,4). The computational basis vectors are related to the eigenvectors & p through 
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where 



eigenvalues are E\ = h(fi mn - a mn ), 



= \Ja 2 mn + SI 2 , A 2 = [pL mn + Q)/4/j, mn , and B 2 = (/j, mn - SY)/4(j, mn . The corresponding 
E 2 = —U (/i mn + o, mn ), E 3 = Ti([i mn + a mn ) , and E4 = h(a mn — /x m „). 



Unitary evolution in truncated basis space 



Unitary time evolution of the system due to interaction operator His had been obtained, analytically, in ref. [8j 
for initial states \g, to — l,n — 1) , and |e, to — 1, n — 1). For interaction time t p such that /z mn i p = pn, p = 1, 2..., the 
initial state \g, m — 1, n — 1) is found to evolve into 



*/(*p) = (-!) P [cos(a mn t p ) \g, m - 1, n - 1) - i sin(a„ m i p ) |e, to, n)] , 



(7) 



Now consider a special initial state with the ion in its ground state occupying the lowest energy trap state, while 
the cavity is prepared in vacuum state (m = n = 1). Without taking into consideration the decoherence effects, the 
ratio a = (nn/an) determines the interaction time t needed to generate maximally entangled tripartite GHZ state. 
Leaving out an overall phase factor, for an interaction time t p such that ant p = ?, (a = 4), the system is 

found to be in the state 



*GHZ,iih) = ^ (Iff. 0, 0> — * |e, 1, 1» , 



and for a\it p = , the state of the system is 



1 



Knzjitp) = ^= (Iff, 0,0) + *|e, 1,1» . 
The density operator representation for these states is 

pGHZ.l{tp) = \*GHZ,l(tp)) (^GHZj{t p )\ . 

Neglecting an overall phase factor, the states can be written in Schrodinger picture as 

*GHz{tp) = ^(lff.O,0)-iexp(-i2wot)|e,l,l)) 
Khz^p) = ^(lfl,0,0) + iexp(-i2wo*)|e,l,l)). 



(8) 



(9) 



(10) 
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DECOHERENCE OF UNITARY EVOLUTION IN TRUNCATED SPACE 

We notice that the interaction Hamiltonian Hn (Eq. connects the model space states to states outside the 
model space. For the specific case of initial state |g, 0,0) , the leading term in the probability of finding a state 
outside the model space varies as t 8 . The states outside the model space may be considered as part of the environment 
states. The resulting decoherence effects can be accounted for by coupling the model space interaction Hamiltonian 
to the environment modeled as a set of oscillators. We recall that the scattering outside the model space results in 
internal states linked to motional states with phonon number > 1, resulting in heating of the ions. In a realistic 
exprerimental set-u p, t he heating of ions trapped in a radio-frequency trap occurs due to the electric field noise from 
the trap electrodes |l7j |. The two effects can not be separated. Error in initial phonon state preparation can also add 
to decoherence by additional scattering to states outside the model space. The cavity states with number of photons 
> 1, are created and annihilated in the vector space outside the model space. In a lossy cavity photons leaving 
the cavity are akin to state reduction which can contribute to creating states in the model space. In this work we 
consider the average decoherence effects arising due to the decoherence sources such as magnetic field fluctuations, 
laser frequency fluctuations, potentials applied to elctrodes and coupling to states outside the model space. The 
random perturbations of the Hamiltonian are accounted for by coupling the model space interaction Hamiltonian to 
environment modeled as a set of non-interacting harmonic oscillators |l2j. The phase decoherence results in decay 
of entanglement needed for various computation related tasks. In interaction picture the Hamiltonian for the system 
and environment, with His coupled to environment is given by 

H = H IS + K"k&t&k + H IS ]T (g* k B k + g k B^j , (11) 

k k 

where B^ and B k are bosonic creation and destruction operators for the environment mode of frequency uj k . The 
coefficients g k , g k are the system environment couplings. 

The initial state of the system is assumed to be ps(0), while the environment is in an uncorrelated state, Jlfc^fc' 
where 9 k = Z k ~ 1 e~P UJk ^k^ k and Z k = (1 — e~^" ,fe ) _1 . The assumption that the environment is in an uncorrelated state 
allows exact solvability, though no physical fundamental reason can be given for this. However, as the initial state of 
the system is a doctored state, the assumption of separability of system environment state at t = 0, holds. The state 
operator for the system is obtained by solving 

Pis(t)=Tr E ^p S (0)n^ f j , (12) 

where Tte refers to the operation of tracing over bosonic degrees of freedom used to model the environment. Working 
in the eigen basis of Hjg, we may write 

i,j \ k / 



PIS 

-i{Ei-Ej)t 



exp 



(Ws(o)\^)\^)(^j\- (is) 



The trace over bosonic bath states is evaluated by using coherent states 0] and the resulting state operator for the 
system is 



p IS (t) = ^(^(O^exp' lE > E ' ]2T[i] 



^3) 

4.0fr 2 



exp 

where, using the notation of ref. Q, 



(Ei - E 3 )t (Ef - E])C{t) 



h 2 



r w- 8 E^f- 2 (T) coth (^)' (15) 
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and 
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C(t) = V i^- (sinKt) - Wfc t) . (16) 



The summation over bath modes in Eqs. (|15|l and Ijltjfl can be replaced by integration over frequency by introducing 



a frequency dependent density function. For an ohmic dissipation characterized by density function D(lo) |g| 2 = 
Kwexp(-w/w c ), we obtain 

r{t) =%K j A /xp^M) co th sin 2 ^fj , (17) 

and 

C(t) = K / da; eX P(-^K) ( sin (^) _ ^) ( 18 ) 

where w c is a cut off frequency and constant k a measure of system-bath coupling strength. 

Initial State ps(0) = \g, m — 1, n — 1} {g, m — 1, n — 1| 

We consider the case when the system is prepared, at t = 0, in the state /3s (0) = |g, m — 1, n — 1) (<?, m — 1, n — 1|. 
The state operator p/s (i) obtained by putting in the energy spectrum of His m Eq. I|14|l is 

PIS® = {A+ 2 B? 0$l)<$ll + l$4) (^1] 

+ (A ~ 5)2 (|$ 3 ) ($ 2 ] + ]$ 3 ) ($3l) 



■ exp 



f ( Mm n n ) ^ ( £ ) ^ 



2 

x [exp (-i2(fi mn - a mn )t) |$i> ($ 4 | + exp {i2((i mn - a mn )t) |$ 4 ) ($i 

(A-B) 2 / . s2 , A 

H ^ exp ^- (/u TO „ + a mn ) r(t)J 



x [exp (i2(fj, mn + a mn )t) |$ 2 ) ($3! + exp (-i2(/x m „ + a m „)i) |$ 3 ) ($ 2 |] 

M 2 - B 2 ) , 2 , „ 
+- g -exp(-/iL r (*)) 

x [exp [-i (2/it m „t - ¥?(*))] |$i) ($ 2 1 + exp [i (2fj, mn t - <p(t))] |$ 2 ) 

- exp [-« (2/w* + I^s) ($4 1 - exp [i (2(i mn t + <p(t))] |$ 4 ) ($3|] 

M 2 - B 2 ) , 2 . 
+- ~ 2 exp (-a 2 mn T(t)) 

x [exp [t (2a mn t + tp(t))] ($ 3 | + exp [-« (2a mn t + p(t))] |$ 3 ) 

- exp [» {2a mn t - <p(t))\ |$ 2 > ($ 4 | - exp [-t (2a mn t - <p(jt))] |* 4 ) ($ 2 |] , (19) 

where y>(t) = 4:/i mn a mn C (t) . The state operator pis{t) contains information about the decay of coherences due to 
coupling of the system with the environment. 

DECOHERENCE AND GHZ STATE GENERATION PROBABILITY 

Undergoing decoherence free time evolution, for ps(0) = \g,0,0) (g, 0, 0| and the interaction time t such that 
dutp = j,(J-iit p = 7r (a = 4), the system evolves into a maximally entangled tripartite two mode GHZ state 

given by Eq. (JSJ. We define GHZ state generation probability as PcHz{t) = tr \Pls(t)p^HZ,i\ where f>&Hz i = 
^ghz i) \ ^ghz i ■ Using Eq. tEH we write down the density operator pis(t) for the choice m = l,n = 1 and 
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evaluate PcHz(t)- The resulting expression for GHZ state generation probability is 

P GH z(t) = l + ^T ( ex P ( "^o^ ) cos(2 m t) cos(y(f)) 

Q 2 ( f-a\ jT(t)\ . ,„ . , n , , 
' exp — 11 v ; sin (2a u t) cos(ip(t)) - 1 



4/ifi V V 4.0 
1 f fin - o u \ 2 / - (/in - an)" 2 



2 Ml i 7 eXP (, 40 jBin[2( m -an)t] 

' ' '" 1 ' " 1 A 2 exp f -(^+«n) 2 r W \ ^ + _ 



2^n / V 4.0 

We also calculate the population inversion defined as I = P g — P e , where P g (P e ) is the probability of finding the ion 
in ground (excited) state. For the choice m = l,n = 1, in the state pis(t) of Eq. lfT§jl. we get 

I = (exp (- (fin - an) 2 T(tfj cos [2(p n - au)t]j 

+ Mii+ffln ( /_ ( + aii)2 r(t) \ cqs [2( + an)t] _ (21) 
2Mn v y 

In the limit k — » the Eqs. (|19|l . (|20|l and i|21|) reduce to decoherence free evolution of the state operator, the GHZ 
state generation probability, and the population inversion, respectively. As the cavity is initially in vacuum state and 
ionic center of mass in the lowest energy trap state, decoherence effects due to cavity decay and heating are expected 
to be small and have not been considered. 



NEGATIVITY AND LINEAR ENTROPY 



The state of three qubit composite system is represented by pis{t) = PABc{t), where A, B, and C refer to the 
internal state of two-level ion, the state of vibrational motion, and the cavity state respectively. A pure state of a 
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tripartite system can be separable, or entangled. Entangled state may or may not have(biseparable) true tripartite 
entanglement. Bipartite entanglement is relatively well understood. Entanglement content of a tripartite system may 
be measured through the entanglement of its bipartite decompositions that is A+(BC), B+(AC), and C+(AB). The 
decoherence effects result in a composite state pABc(t), which is usually a mixed state. Mixed state entanglement is 
less understood than the pure state entanglement. On examining Eq. I|19fl . m which the state operator is expressed in 
the eigen-basis of His, we notice that the diagonal matrix elements of density matrix are not affected by coupling to 
the environment, whereas the moduli of off diagonal matrix elements decay with time. Consequently, when only non- 
dissipative decoherence is present the trace of PABc{t) remains constant with time and the ei genv alues continue to be 
positive. Since pABc{t) is positive definite, we use negativity proposed by Vidal and Werner as an entanglement 
measure for pure as well as mixed states. Negativity is based on partial transpose of density matrix of composite 
system with respect to a subsystem. 

The state of a bipartite system, composed of subsystems A and B in finite dimensional Hubert spaces cLa and cIb, 
is written as 

dA dn 

p= ( i ' m \p\3,r)\ i,m) (j,r\. (22) 

ij — l m,r— 1 

The partial transpose of density operator with respect to sub-system A is defined as 

d A d B 

P TA = E E < z < m Id 3, r) I j, m) (i,r\. (23) 

ra,r 

The partial transpose of density matrix of an entangled state is not positive definite. Negativity defined as N A — 
y~) j | A, | , where are the negative eigenvalues of p TA is a measure of entanglement of quantum system A with 
B. For a separable state N A — and at maximal entanglement value of N A depends on dimension dA- For the 
tripartite system at hand, we construct the partial transpose with respect to sub-systems A, B or C, while keeping 
the composite state of the remaining two sub-systems unaltered. For the density matrix operator pABc(t) acting on 
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k=0.00, k=0.001 

k=0.01, k=0.02 

K =0.05, -k=0.1 




FIG. 3: N, calculated from p TA , as a function of scaled time variable T(— ant), for the choice /in\an = 4, and k = 
0, 0.001, 0.01, 0.02, 0.05 and 0.1. 

composite space, the transpose with respect to sub-system A reads as 

dA ds dc 

P Ta = $Z ^2 (i,m,n\p\ j,r,s)\ j,m,n) (i,r,s\, (24) 

— l m.r—1 n,s — 1 

where dx{X = A,B,C) refers to dimension of the Hilbert space of subsystem X. We use p TA , p TB , and p Tc to 
calculate the negativitiy measures for sub-systems A, B, and C , respectively 

For a quantum state p in d dimensional Hilbert space, the linear entropy Si is defined as 

S, = ^4i(l-Tr(^)). (25) 

Linear entropy is used to measure the mixedness of subsystems A,B, and C. For pure states Si = 0, whereas 
for maximally mixed states Si = 1. Reduced state operators p^ ed = trsciPABc), Pred = t r Ac(PABc), and p^ ed 
= trAsiPABc) are used to calculate linear entropy Si, for the subsystems A, B, and C, respectively. Negativities 
and linear entropies are used to characterize the state of the composite system at an instant t, for a given value of 
coupling to the environment. 

NUMERICAL RESULTS AND CONCLUSIONS 

We consider the system prepared in the state /5s (0) = \g, 0, 0) (g, 0, 0|, at t = 0. The laser-ion coupling constant is 
taken to be f2 = 8.95 MHz [5j, where as the ratio fin /an — 4.0. Numerical values of T(t) and C(t) are obtained 
by solving the integrals of Eqs. (|17I18|I for the choice uj c = 1200MHz, temperature T = 0.03K and different values 
of coupling strength k. The variable w in Eqs. (|17I18|I is varied from zero to 3w c . The value k = 0.001 corresponds 
to a weak coupling and successively larger values tend to produce decoherence on a scale comparable to GHZ state 
generation time which is of the order of 0.34 nano s. The state of the system at an instant t can be detected 
experimentally by cavity- photon measurement combined with atomic population inversion measurement. 

Figs. (1) and (2) display PaHz{t) and population inversion I(t), respectively, as a function of scaled time variable 
T(= ant) for decoherence parameter values of k = 0,0.001,0.01, and 0.1. Evidently the evolution dynamics of 
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FIG. 4: N, calculated from p 5 ^ = p Tc , as a function of scaled time variable T(= ant), for the choice pn\an = 4, and 
k = 0, 0.001, 0.01, 0.02, 0.05 and 0.1. 

PoHzf) and If) is sensitive to changes in the decoherence parameter k. The presence of decay factors in Eq. I)20|) 
causes the peak value of Pghz to decrease with time. Phase decoherence causes a shift in the interaction time 
needed to get the system in maximally entangled state. The population inversion however does not show any phase 
decoherence effects. We may point out that at T = 135° the tripartite system is in maximally entangled state of Eq. 
@ which is orthogonal to the state of Eq. (JSJfor which Pghz if) has been calculated. Population inversion is zero at 
T = 45° as well as 135°, signalling maximally entangled state generation. For weak coupling (k = 0.001), the behavior 
of If) does not differ much from no decoherence curve. However, for larger values of n, If) for the separable system 
at T = 45° can indicate the amount of decoherence present. 

The state operator of Eq. i|19[l for m = 1, n = 1 is used to calculate the density matrices transposed with respect 
to quantum systems A, B, and C. The negativities are calculated, numerically, from p TA , p TB , and /? Tc and plotted as 
a function of scaled time variable T{= ant). The negativity plot displaying entanglement dynamics of internal state 
of two-level ion (A) is shown in Fig. 3. For subsystem B, the negativity calculated from p TB =p Tc , is displayed as a 
function T(= ant) in Fig. 4. The values of decoherence parameter are k = 0, 0.001, 0.01, 0.02, 0.05 and 0.1. The time 
evolution of vibrational motion of ion's center of mass(_B) and cavity number state (C) are identical due to inherent 
symmetry of these subsystems. We notice that as the value of decoherence parameter increases, there is a decrease in 
the maximum entanglement of sub-systems A, B and C with their compliments. Besides that the phase decoherence 
changes the time at which maximal entanglement of the three parties is observed. These two factors are consistent 
with the behavior of Panzf) and If) in Figs. 1 and 2. 

Linear entropies Si, for subsystems A, B and C have been calculated using the reduced density operators obtained 
from pisf) of Eq. i|19|) f° r m — \,n = \. Figs. 5 and 6, display Si as a function of scaled time variable T(= ant), 
for subsystems A and B(or C), respectively. The decoherence parameter values are k = 0,0.001,0.01,0.02,0.05 and 
0.1. For k = 0.001,0.01, the trend of time evolution of linear entropy follows the variation of entanglement as seen 
in negativity plots of Figs. 3 and 4. However, for a stronger coupling value of k = 0.02 for subsystem A, whereas JV 
tends to zero for T approaching 90°, Si shows fluctuation for T > 90°. Similarly for k — 0.05 negativity indicates 
a decoupling of system A from system BC beyond T = 67° while Si points to the presence of some correlations. 
For k = 0.1, the system is highly damped and both JV and Si show the ion to become separable around T = 45°. 
We may conclude that for the states lying at the boundary of entangled separable region the negativity and linear 
entropy show distinctly different variation with time the reason being that the negativity measures only entanglement 
generating quantum correlations, whereas the linear entropy measures all correlations that reduce the purity of the 
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state. 

To conclude, we have investigated the decoherence process of three qubit system obtained by manipulating the 
state of a cold trapped two-level ion coupled to an optical cavity. Interaction of the ion with a resonant laser field 
and the cavity field tuned to red sideband of ionic vibrational motion, generates entanglement of internal states 
of the ion, vibrational states of ionic center of mass and the cavity field state. Non-dissipative decoherence of the 
state of the system occurs due to interaction of the system with the environment. We have considered the effect of 
decoherence sources such as magnetic field fluctuations, laser frequency fluctuations, potentials applied to elctrodes 
and coupling to states outside the model space. The random perturbations of the Hamiltonian are accounted for 
by coupling the model space interaction Hamiltonian to the the environment, modeled as a set of non- interacting 
harmonic oscillators. The pointer observable is energy of the isolated quantum system. Analytic expression for the 
state operator of tripartite composite system, including non dissipative decoherence effects has been obtained. The 
initial state of the tripartite system is a separable state and the state of environment has no initial correlations. 
Coupling to environment is found to introduce an exponential decay with time of the off diagonal matrix elements 
of density operator pis(t). In addition, a defasing of various states in the superposition also occurs. The extent to 
which the process is decohered depends on the system envionment coupling strenth, that can be partially controlled 
by adjusting the electrode potentials. Besides the strength of coupling to the environment, the energy spectrum 
of system Hamiltonian plays an important role in determining the decoherence rate of a given initial state of the 
composite system. For a specific choice of interaction parameters, the isolated system evolves to tripartite GHZ 
state 8]. We have evaluated analytically the probability of generating maximally entangled GHZ state, and the 
population inversion in the presence of non-dissipative decoherence. Numerical calculations for different values of 
system environment coupling strengths using interaction parameter from a recent experiment show the peak value 
of GHZ state generation probability to decrease with increase in k. PcHz(t) is also sensitive to phase decoherence. 
The population inversion however does not show any phase decoherence effects. The results can be used to study 
the effects of engineered environment on the process of tripartite maximally entangled state generation. Bipartite 
entanglement of cavity mode and atomic internal states can be produced without coupling to the motional states. The 
tripartite entanglement generated through coupling to motional states can be transferred to other ions that might be 
added to the trap. In comparison with an earlier proposal for GHZ state generation |13|, in the current scheme the 
initial state of the system is \g,0,0). As such the cavity losses and spontaneous decay effects are minimized. 

Using negativity as an entanglement measure, the entanglement dynamics of the tripartite system in the presence 
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FIG. 6: Linear entropy Si, calculated from pf etj , as a function of scaled time variable T(= ant), for the choice /ttn\on = 4, 
and k = 0, 0.001, 0.01, 0.02, 0.05 and 0.1. 

of decoherence has been analyzed. As the value of decoherence parameter increases, the maximum entanglement of 
sub-systems A and B with their compliments decreases. In the context of trapped ion radiated by the single mode 
cavity field and an external resonant laser, the time for which the two level ion remains coupled to the state of 
vibrational motion and the cavity state decreases with increase in k. For large values of k the states of the composite 
system lying at the boundary of entangled-separable region are reached. For these states the negativity and linear 
entropy show distinctly different variation with time. This can be understood by noting that whereas the negativity 
measures entanglement generating quantum correlations, the linear entropy measures all correlations that reduce the 
purity of the state. 
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